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1  Introduction 


■  This  report  will  address  the  effect  of  a  ‘^bump’’^in  the  turbulence  spectrum  on 
the  propagation  of  a  laser  beam.  The  motivation  for  this  work  is  as  follows:  In  con¬ 
vent!  mal  treatments  on  the  effect  of  turbulence  on  laser  propagation,  it  is  usually 
assumed  that  the  turbulence  can  be  described  by  a  spectral  density  which  represents 
a  cascade  of  energy  from  long  wavelength  modes  down  to  shorter  wavelength  modes. 

The  most  common  assumption  is  that  of  the  Kolmogorov  models  where  in  the  ‘'in¬ 
ertial”  range  of  the  spectrum  (between  a  minimum  wavenumber  and  maximum 
wavenumber  At/q,  corresponding  to  the  maximum  size  Lq  and  minimum  size  Iq  of  the 
turbulent  eddies)  the  spectral  density  decreases  as  ^  dimensions).  In 

some  media  however,  a  source  for  injecting  energy  into  the  turbulent  spectrum  at  a  ,  ^ 

particular  discrete  wavelength  (or  over  a  range  of  wavelengths)  may  be  present:  for  ' 
example,  this  is  the  case  in  a  turbulent  plasma  where  an  instability  may  be  grow¬ 
ing  over  a  range  of  wavelengths.  The  presence  of  such  an  energy  source  will  thus 
produce  a  “bump”  in  the  turbulence  spectrum^  as  shown  schematically  in  Fig.l.  As 
we  shall  see  below,  the  spectrum  can  then  be  viewed  as  this  “bump”^superimposed 
upon  the  background  natural  cascade;  that  is,  the  spectrum  can  be  analyzed  as  the 
sum  of  a  cascade  spectrum  and  the  “bump”^spectrum,  and  the  effects  of  each  on  the 
propagation  of  the  beam  can  thus  be  analyzed  separately.  (  C  '  ^  — - 

This  report  is  organized  as  follows.  In  Section  2  we  will  briefly  review  the  math¬ 
ematical  formalism  that  will  be  used  to  analyze  the  statistical  features  of  the  effect 
of  turbulence  on  laser  propagation.  In  Section  3,  this  formalism  will  be  applied  to 
the  case  where  the  turbulence  is  described  by  a  spectral  density  which  is  a  “bump” 
at  some  finite  wave  number  kb.  Section  3.2  will  treat  the  idealized  (but  analytically 
tractable)  case  where  the  “bump”  is  a  5-function;  the  effect  of  a  “bump”  of  finite 
width  is  studied  numerically  in  Section  3.3.  Finally,  the  superposition  of  a  “bump” 
on  a  cascade-type  spectrum  is  considered  in  Section  3.4.  One  important  conclusion 
which  results  from  this  study  is  that  in  three  dimensions,  even  if  the  amplitude  of 
the  “bump”  component  of  the  spectrum  is  much  smaller  than  the  amplitude  of  the 
cascade  component  (for  example,  by  ten  orders  of  magnitude),  if  the  scalelength  of 
the  turbulence  represented  by  the  “bump”  is  much  smaller  than  the  scalelength  of 
the  cascade  turbulence  (say  four  orders  of  magnitude  smaller),  then  almost  all  f  the 
power  in  the  spectrum  is  due  to  the  “bump”  component.  Furthermore,  in  this  case, 
the  effects  of  the  turbulence  on  the  beam  are  in  effect  determined  by  the  “bump”  com¬ 
ponent,  and  the  existence  of  the  ccLscade  component  (which  may  represent  turbulence 
at  very  large  scalelengths)  can  be  ignored. 
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Figure  1.  Example  of  a  turbulence  spectrum  with  a  “bump”  superimposed 
on  a  cascade. 
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2  Preliminaries 


In  this  Section  we  will  briefly  review  the  mathematical  constructions  required  to 
analyze  the  statistical  effects  of  a  turbulent  medium  on  the  propagation  of  a  laser 
beam.  For  simplicity,  we  will  assume  that  the  problem  is  two-dimensional:  the  beam 
(with  wavenumber  ko)  is  taken  to  be  propagating  in  the  ^-direction,  with  a  transverse 
beam  profile  ip{x)  on  the  i-axis.  The  extension  of  mathematical  expressions  to  three 
dimensions  is  straightforward,  and  the  main  results  are  not  changed.  Furthermore, 
we  shall  assume  throughout  this  article  that  the  laser  wavelength  A  =  2Tr/ko  is  much 
shorter  than  any  scalelength  of  variation  of  the  medium  (f.e.,  the  minimum  size  lo  of 
the  turbulent  eddies)  so  that  we  are  working  in  the  eikonal  (or  WKB)  regime. 

If  the  disturbance  of  the  beam  due  to  the  turbulent  medium  is  small  over  a  propa¬ 
gation  distance  z,  the  beam  can  be  thought  of  as  a  bundle  of  rays  propagating  roughly 
parallel  to  the  2-direction.  The  phase  shift  <j)  along  a  ray  propagating  at  the  trans¬ 
verse  location  x  caused  by  the  turbulent  perturbations  8N{x^z)  in  the  media’s  index 
of  refraction  is  therefore  ^ 

(f>{x)  =  ko  f  dz'  SN(x,z')  (1) 

Jo 

In  this  expression,  4>  is  the  phase  shift  induced  by  a  single  realization  of  the  turbulence, 
characterized  by  6N. 

2.1  Turbulence  and  Phase  Correlation  Functions 

A  statistical  analysis  of  the  effect  of  the  turbulence  on  beam  propagation  is  based 
on  averages  of  characteristic  quantities  (such  as  this  phase  shift)  over  an  ensemble 
of  realizations  of  the  turbulence.  Denoting  this  average  by  the  angle  brackets  (•), 
we  cLssume  that  the  turbulent  fluctuations  have  zero  mean,  {6N{x,z))  =  0,  so  that 
the  mean  phase  shift  from  (1)  also  vanishes.  The  correlation  function  of  the  phase  is 
defined  as 

R^{\x  -  x'l ;  z)  =  {<P{x)<f>{x'))  =  kl  r  dz'dz"  {SN[x,z')8N{x\z")) 

Jo 

=  kl  J  dz'  dz"  Ri\f  ^^/(x'  —  xy  +  {z"  —  z')'^  (2) 

Here  we  have  assumed  that  the  turbulence  is  isotropic  and  homogeneous:  the  re¬ 
fractive  index  correlation  function  iZyv  depends  only  the  the  spatial  separation  of  the 
two  points  (x,2')  and  (x',  2").  In  this  case,  the  phase  correlation  function  depends 
only  on  the  separation  |x  —  x'\  on  the  transverse  x-axis  (and  the  propagation  distance 
z).  Because  of  this  assumption,  this  expression  can  be  simplified:  changing  variables 
Sx  =  x'  —  X,  =  2"  —  2',  it  can  be  recast  as 

RM  ;  =  kl  £  dz'  ds,  Rr,  +  52)  (3) 
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The  area  in  (2',  52)-space  over  which  the  integral  is  to  be  performed  is  shown  in 
Fig. 2a;  here,  the  integral  in  the  s^-direction  is  to  be  performed  with  z'  fixed,  and 
then  z'  is  varied.  Changing  this  order  of  integration  requires  changing  the  limits  of 
integration:  holding  Sj  fixed,  first  integrate  z'  from  — to  2  (for  <  0),  or  from  0 
to  2  —  (for  Sj  >  0),  and  then  vary  from  —2  to  2  (see  Fig.2b). 

^0(1-Sx|;2)  =  ^0  j  j  dz'  Rn  {\[sl  +  52^ 

+kl  £  ds,  dz'  Rn  (\AF+^)  (4) 

The  integrals  over  z'  can  now  be  performed;  then  letting  — >  —s^  in  the  first  integral, 

the  two  integrals  can  be  combined  to  give 

R^{\sx\]z)  =  2kl  dsz  (2  -  s^)  Rn  (^\/sl  +  (5) 

We  note  that  this  is  an  exact  expression  for  the  relationship  between  the  correla¬ 
tion  function  R^  of  the  wave  phase  and  the  correlation  function  Rn  of  the  refractive 
index  fluctuations,  ba^ed  only  on  the  definition  (1)  for  the  phase  shift.  In  practice, 
however,  Rn  typically  decays  to  zero  over  some  correlation  distance  Lo  (often  called 
the  outer  scale  of  the  turbulence),  and  this  distance  is  much  less  than  the  propa¬ 
gation  distance  of  interest  {Lq  <C  ^).  In  this  case,  the  maximum  value  of  that 
contributes  to  the  integral  (for  fixed  Si,)  is  about  Lo,  so  that  2  —  «  2  is  a  very  good 

approximation.  Taking  this  term  outside  the  integral  then,  we  have 

/?,^(|sxl ;  z)  «  2kQZ  dsz  Rn  (6) 

where  the  same  approximation  allows  us  to  extend  the  limit  of  integration  to  infinity. 

2.2  Turbulence  and  Phase  Spectral  Densities 

The  spectral  density  of  a  homogeneous,  isotropic  random  process  (such  as  the 
turbulent  index  of  refraction  or  the  induced  phase  shift  of  the  Iciser  beam)  is  defined  as 
the  Fourier  transform  of  the  correlation  function.  For  the  two-dimensional  turbulence 
we  have 

SN{k,,k,)  =  J  ds^ds,  Rn  [y/si+sl^ 

sds  Rn{s)  / 

Jo 

SN{k)  =  2Tr  J  sds  Rn{s)  Jo{ks) 


(7) 
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Figure  2.  Area  over  which  (2',  integrals  to  be  performed  in  (3) 
(a)  Original  order  of  integration,  (b).  Inverted  order  of  integration 
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Since  the  turbulence  is  homogeneous  and  isotropic  in  physical  space,  it  is  isotropic  in 
/c-space  as  well.  The  inverse  of  (7)  is  similar; 


In  practice,  the  spectral  density  is  often  the  more  fundamental  quantity  in  the 
sense  that  a  turbulence  model  is  typically  specified  by  postulating  the  form  of  the 
spectral  density  5/v-  One  reason  for  this  is  that  it  can  be  shown  that  the  spectral 
density  of  a  physically  realizable  random  process  must  be  non-negative.  Thus,  spec¬ 
ifying  a  model  for  the  correlation  function  /2/v  may  produce  by  construction  (7)  a 
spectral  density  which  does  not  satisfy  this  criterion;  beginning  with  a  form  for  Sn 
avoids  this  problem,  and  the  correlation  function  can  then  be  constructed  from  (8). 

The  spectral  density  associated  with  the  one-dimensional  phase  correlation  func¬ 
tion  is  defined  by 

St(k.)  =  f  ds,  e-''-" 

RtM  =  /  ^  e*"  (9) 

Note  that  since  the  phase  fluctuations  are  a  one-dimensional  random  process  (in 
this  two-dimensional  propagation  model),  the  Fourier  transforms  in  (9)  are  one¬ 
dimensional  (as  opposed  to  those  in  (7,8)  for  the  two-dimensional  random  process 
of  the  refractive  index  fluctuations). 

Since  the  phase  correlation  function  and  refractive  index  correlation  function  are 
related  by  (5),  a  relation  between  their  a.ssociated  spectral  densities  can  ?.iso  be  found. 
This  is  most  easily  achieved  by  introducing  the  first  line  of  (8)  into  (2)  to  obtain 

R^{S:r)  =  kl  dz'  dz"  (\/sl  +  (z"  - 

=  ki  r  dz'dz"  f  (10) 

Jo  J  ZTi  Zw 
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Finally,  performing  the  Fourier  transform  of  both  sides  to  obtain  S^{kx)  (as  in  (9)) 
the  result  is 

=  (5^)'  (12) 

This  expression  could  be  obtained  beginning  with  (5)  as  well. 

The  interpretation  of  (12)  is  that  the  spectral  density  of  the  phase  fluctuations 
on  the  transverse  x-axis  is  just  the  spectral  density  of  the  two-dimensional  refractive 
index  fluctuations  projected  onto  the  x-axis,  but  filtered  so  that  only  turbulent  pertur¬ 
bations  with  wavelengths  in  the  direction  of  propagation  (the  z-direction)  comparable 
to  the  propagation  distance  z  contribute.  Two  observations  are  immediate.  First,  if 
the  spectral  density  of  the  turbulence  is  non-negative  (as  it  must  be  to  describe 
a  physically  realizable  random  proce.ss),  then  (12)  clearly  shows  that  the  phase  spec¬ 
trum  S4,  is  non-negative  (note  that  this  relation  follows  exactly  from  the  premise  (1) 
for  the  wave  phase  shift).  This  fact  is  extremely  important  in  practice,  where  the 
so-called  phase-screen  method  is  often  used  to  numerically  propagate  the  laser  beam 
through  a  turbulent  medium.  In  that  method,  the  wave  is  propagated  as  though  in  a 
uniform  non-turbulent  medium  a  distance  z  between  “phase-shifting  screens”  which 
represent  the  effect  of  the  turbulence;  on  each  “screen”,  the  wave  V>(x)  is  modulated 
by  a  factor  exp(i<^(x;  z)),  where  </»(x;  z)  is  a  random  phase  function.  The  phase  func¬ 
tion  <t>{x]  z)  is  generated  by  first  constructing  its  Fourier  transform  z)  from  a 
Gaussian  random  distribution  with  variance  S^{kx)  at  each  kx,  and  then  transform¬ 
ing  this  into  x-space.  It  is  straightforward  to  show  that  this  method  will  produce  the 
correct  phase  statistics  (i.e.,  the  specified  spectral  density  and  correlation  function 
of  the  phase)  over  an  ensemble  of  a  large  number  of  phase  screens.  This  method 
therefore  inherently  requires  that  S<i,{kx)  be  non-negative,  since  its  value  at  each  kx  is 
used  as  the  variance  of  a  probability  distribution.  As  we  have  seen  above,  this  qual¬ 
ity  is  guaranteed  by  first  specifying  the  spectral  density  Sti  of  the  refractive  index 
fluctuations  to  be  non-negative,  and  then  constructing  from  (12). 

The  integral  in  (12)  presents  no  particular  problem  for  numerical  computation, 
and  in  fact  should  be  quite  easy  to  evaluate  in  most  cases  of  interest  since  the  “filter 
function”  [sin  A:jz/2/(A;jz/2)]^  falls  off  rapidly  (its  first  zero  is  at  k^  =  2ttIz).  Our 
second  observation,  however,  is  that  if  5//  is  slowly  varying  in  the  vicinity  of  =  0 
(for  all  kx)  then  (12)  can  be  approximated  by  setting  =  0  in  the  argument  of  S'yv 
and  performing  the  integral  only  over  the  “filter  function”  to  obtain 


S4,{kx)  =  klz  S^ikx.kx  =  0)  (13) 

Again  we  see  that  5,^  >  0  if  Sjw  is.  The  only  requirement  for  (13)  to  hold  is  that 
varies  slowly  in  the  A:^-direction  in  the  vicinity  kzzl2  <  tt  for  all  values  of  kx\ 
thus,  only  wavelengths  in  the  direction  of  propagation  comparable  to  z  are  admitted, 
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and  their  contribution  is  approximated  by  a  flat  distribution  at  the  —  0  value.  It 
is  furthermore  quite  simple  to  show  that  while  the  inverse  Fourier  transform  of  the 
exact  relation  (12)  yields  (5)  (as  expected),  the  inverse  transform  of  the  approximate 
relation  (13)  produces  the  approximation  (6);  thus,  the  approximations  leading  to 
both  (6)  and  (13)  from  (5)  and  (12),  respectively,  are  consistent. 

As  a  final  note,  one  can  relate  the  level  of  the  phase  fluctuations  to  the  level  of 
the  refractive  index  fluctuations.  These  quantities  (which  have  been  .assumed  to  be 
independent  of  location)  are  defined  as 

{8N^)  =  i?yv(0)  (14) 

(<^2)  =  R40)  (15) 

Using  (8)  and  (9),  the  value  of  either  correlation  function  Rj\j  or  R^  at  the  origin  j 
given  by  an  integral  over  the  associated  spectral  density.  Thus,  using  the  approximate 
form  (13)  we  have 

f  (Ik 

=  Rm  =  J^  s,{k^) 

=  klz  J  k,  -  0) 

=  klz{8N^)ii  (16) 

The  one-dimensional  integral  over  defines  the  integral  scalelength  £{  of  the  turbu¬ 
lence. 


2.3  Average  Intensity 

We  conclude  this  Section  with  a  discussion  of  the  effect  of  a  random  medium 
(characterized  by  a  spectral  density  for  the  refractive  index  fluctuations)  on  the  in¬ 
tensity  profile  of  a  laser  beam.  In  the  absence  of  turbulence  (or  other  effects  of  the 
medium,  such  as  thermal  blooming),  a  coherent  beam  of  initial  radius  Tq  will  spread 
due  to  natural  diffraction  according  to 


^0 


(17) 


where  Lr  =  kor^  is  called  the  Rayleigh  length  of  the  beam.'  This  spreading  can  be 
viewed  as  the  coherent  interference  of  a  bundle  of  rays  propagating  within  en  angle 

'This  expression  is  for  two-dimensional  propagation  (in  three  dimensions,  the  factor  of  1/4  be¬ 
comes  unity).  In  this  way,  Lr  is  the  same  in  both  two  and  three  dimensions,  for  a  given  initial  radius 
(or  half-width)  tq  (defined  by  (26)  in  both  cases);  the  difference  then  is  that  in  two  dimensions  the 
width  will  diouble  in  a  distance  c  =  y/vlLR,  while  in  three  dimensions  the  area  will  double  at  i  =  Lr. 
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0  —  k±fko  ~  ('toTo)"'  of  the  2r-axis,  where  the  transverse  wavenumber  k±  is  inversely 
proportional  to  the  initial  beam  radius.  In  a  turbulent  medium,  however,  we  see  from 
(1)  that  the  phase  along  each  ray  develops  a  random  component  depending  on  its 
path  through  the  randomly  fluctuating  index  of  refraction.  Thus,  if  the  level  of  the 
refractive  index  fluctuations  is  large  enough,  or  the  propagation  distance  long  enough 
(so  that  the  phase  fluctuations  are  large  enough),  the  effects  of  these  random  phase 
shifts  in  the  the  interference  of  the  rays  will  begin  to  be  seen.  Two  major  effects  are  a 
loss  of  transv'erse  beam  coherence,  and  beam  broadening  due  to  the  turbulence.  The 
loss  of  coherence  is  manifested  by  a  loss  of  smoothness  in  the  transverse  intensity 
profile;  the  intensity  develops  “spikes”  so  that  it  appears  to  have  broken  up  into 
smaller-diameter  “b“amlets”.  Beam  broadening  can  be  thought  of  either  as  the  result 
of  a  non-coherent,  mostly  forward  scattering  of  the  rays  off  of  the  turbulent  refractive 
index  fluctuations,  or  as  a  result  of  the  more  rapid  natural  broadening  of  the  smaller- 
diameter  beamlets  (with  smaller  Rayleigh  lengths)  produced  by  the  loss  of  coherence. 

Due  to  the  randomness  in  the  turbulence,  it  is  impossible  to  predict  the  evo¬ 
lution  of  the  laser  intensity  profile  as  it  propagates  through  a  single  realization  of 
the  turbulence.  However,  over  many  realizations  of  the  turbulence  governed  by  a 
spectral  density  (or,  over  many  laser  shots  through  an  ever- fluctuating  medium), 
predictions  can  be  made  regarding  the  average  behavior  of  the  beam  intensity  and  its 
characteristics  (such  as  its  width  or  coherence  length).  For  example,  a  prediction  for 
the  ensemble- averaged  wave  intensity  {I{x,z))  =  {\rp{x,z)f)  can  be  derived^  based 
on  a  modified  Huygens  principle: 

{I{x,z))  I  ds  MT{sl2;z)MA{s;z)e-'^°^^^^  (18) 

Here,  the  turbulence  modulation  correlation  function  Mx{s]  z)  is  defined  as 

Mt{s;  z)  =  (19) 

while  the  amplitude  correlation  Ma{3]z)  is 

Ma{s-,z)  =  j  dx  xl>{x  +  \s,t))xf{x-\s,t))e''^°=^^l^  (20) 

It  can  be  shown'  that  if  the  random  phase  shift  (f){x,z)  is  a  Gaussian  random 

variable  at  each  x,  then  the  turbulence  modulation  correlation  function  becomes 

A/7’(s;  z)  =  (21) 

where  (0^)  is  the  level  of  the  phase  fluctuations  (16)  and  R:f,{s)  is  the  phase  correlation 
function.  We  note  that  since  Mt(0;  2)  =  1,  the  total  power  in  the  beam  {i.e.,  the 
integral  of  (18)  over  .r)  remains  constant  as  it  propagates  in  z. 
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For  an  initial  Gaussian  amplitude  profile 


=  (22) 

the  amplitude  function  Ma{s]z)  is 

Ma{s-,  z)  =  (23) 

where 

Note  that  in  the  absence  of  turbulence  ((<^^)  =  0),  Mt{s;z)  is  unity  so  that  (23)  in 
(18)  yields  the  correct  expression  for  the  propagation  of  a  Gaussian  pulse  in  vacuum: 


I{x,z) 


_Z_ 

a{z) 


(25) 


Furthermore,  (24)  clearly  shows  the  effect  of  natural  diffraction  as  stated  in  (17),  with 
r(z)  =  <t(2)/v^. 

One  use  of  the  expression  (18)  is  to  compute  the  degree  of  beam  broadening  due 
to  turbulence.  Defining  the  width  r  of  the  beam  to  be 


r\z) 


J  dx  x^  {I{x]  z)) 

j  dx  {I{x\z)) 


(26) 


one  can  show^  that  the  width  grows  as^ 

r'^{z) 


rl 


Here,  the  turbulent  broadening  length  Lt  is  defined  as 

I  1/3 


Lt  = 


rie, 


(27) 


(28) 


in  terms  of  the  level  of  refractive  index  fluctuations  {SN"^),  the  integral  scalelength  £i 
of  the  turbulence  (16),  and  the  transverse  scalelength  of  the  turbulence  £t'- 


J  dk  k'^Sr^ik) 

j  dk  SN{k) 


(29) 


f  As  in  (17),  the  factor  of  1/4  becomes  unity  in  three  dimensions,  and  the  length  Lt  is  the  same 
in  both  two  and  three  dimensions  for  a  given  radius  (or  half-width)  tq. 
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Comparing  (27)  with  (17),  we  see  that  if  the  turbulent  broadening  length  is  compa¬ 
rable  to  or  less  than  the  Rayleigh  length,  then  turbulent  broadening  will  dominate 
that  due  to  natural  diffraction,  and  the  beam  will  broaden  much  more  rapidly  with 
the  propagation  distance. 

The  relation  (18)  describes  the  mean  evolution  of  the  laser  intensity  (averaged  over 

many  laser  shots  through  different  realizations  of  the  turbulence).  An  expression  can 

2 

also  be  derived^  for  the  mean  evolution  of  the  wave  spectrum,  {I{k,  z))  =  (  z)  ): 

{I{k,  ^  0)1'  Mrik  -  z)  (30) 


where  iplk^z)  is  the  Fourier  transform  of  ‘>p{x,z)  and  My  is  the  Fourier  transform  of 
My 

My{k-  z)  =  I  ds  My{s;  z)  e-’''*  (31) 


In  the  absence  of  turbulence,  My(s;z)  =  1  so  that  My{k\z)  =  27r6(fc);  thus  we 
see  that,  a.s  expected,  the  wave  spectrum  from  (30)  does  not  change  as  the  wave 
propagates  in  vacuum.  Scattering  of  the  wave  by  turbulent  fluctuations  will  produce 
changes  in  the  spectrum,  primarily  by  causing  a  shift  to  larger  and  larger  transverse 
wavenumbers;  this  then  is  the  mechanism  for  a  broadening  of  the  wave  spectrum  due 
to  turbulence. 


3  The  Effect  of  a  “Bump”  in  the  Spectrum 

In  this  Section  we  will  use  the  mathematical  formalism  of  the  preceding  Section 
to  analyze  the  effect  of  a  “bump”  in  the  turbulence  spectrum  on  the  propagation  of 
a  laser  beam. 


3.1  Discussion  of  Superposition 

We  first  assume  that  we  are  given  an  isotropic  two-dimensional  turbulence  spec¬ 
trum  S!\f{k)  (depending  only  on  the  magnitude  of  k  =  \jkl  +  k'^)  similar  to  that 
shown  in  Fig.l.  Such  a  spectrum  can  be  assumed  to  be  a  superposition  of  a  “bump” 
spectrum  5^^  on  a  natural  cascade  spectrum  S^\  or 

5yv(fc)  =  4  (32) 

Since  all  the  expressions  involving  the  spectral  density  in  the  preceding  Section  (such 
as  the  Fourier  transforms  in  (7-9)  and  the  relations  between  the  phase  and  turbulence 
spectra  (12,13))  are  linear  in  the  spectral  density,  it  would  appear  that  we  could  treat 
each  component  of  the  spectrum  (5^^  and  5^^)  independently  and  then  linearly 
superpose  the  results. 

Consider,  however,  what  this  means  in  terms  of  the  numerical  phase-screen  method. 
As  briefly  described  in  Section  2.2,  this  method  is  based  on  first  constructing  a  random 
phase  shift  4>{k)  at  each  value  of  k  from  a  Gaussian  random  probability  distribution 
with  the  value  of  S^{k)  as  the  variance.  The  justification  for  this  comes  from  the 
following  considerations.  From  (9)  we  have  (with  kx—*k) 

-  x'l)  =  (*{x)«(i'))  =  /  ^  s*(t)  (33) 

Now  performing  the  Fourier  transform  over  x  and  x'  we  obtain 

{^{k')^‘{k"))  =  I  dxdx'  J  ~  S^{k) 

=  2nS4k')6{k' -  k")  (34) 

This  shows  that  a  homogeneous  random  process  (such  as  the  phase  (f>{x))  is  charac¬ 
terized  by  a  Fourier  transform  <f){k)  which  is  an  independent  random  variable  with 
variance  S4,{k)  at  each  k. 

Now  suppose  we  have  the  case  as  in  (32)  where  the  spectrum  (constructed  from 
S/v  by  either  of  the  linear  relations  (12)  or  (13))  is 

S*(t)  =  Si'4)  +  (35) 
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To  consider  the  effects  of  each  component  of  the  spectrum  separately  and  then  su¬ 
perpose  the  results  means  that  we  should  be  able  to  separately  construct  a  phase 
shift  (i>c{x)  (from  using  S^\k)  as  a  variance)  and  a  phase  shift  (f>bix)  (from 

using  S^^\k)  as  a  variance)  and  then  add  them  together  for  a  total  phase  shift 
(f){x)  =  (i)c{x)  +  (t>b{x)  which  has  the  statistics  (over  an  ensemble  of  phase  screen  real¬ 
izations)  of  the  total  spectrum  S^.  In  A:-space  this  means  that  (34)  must  hold  with 
^{k)  =  (i>c{k)  -f  (pbik)  and  given  by  (35).  Expanding  the  left-hand  side  gives 

{m^-{k'))  =  mk)  +  Mk)U:ik')  +  m')]) 

=  {Uk)m'))  +  {k{k)k{k')) 

Hkmi{k‘))  +  {k{k)k{k')) 

+  iMkmk'))  (36) 

Thus  we  see  that  for  the  effects  of  each  component  of  the  spectrum  to  be  treated 
separately  and  then  superimposed  requires  that  the  cross-correlation  of  the  phase 
shifts  generated  by  each  component  of  the  spectrum  vanishes.  This  means  that  we 
must  assume  much  more  than  simply  that  the  spectrum  can  be  decomposed  into  the 
sum  of  two  (or  more)  components  as  in  (32)  or  (35)  (as  any  spectrum  obviously  can 
be);  we  must  furthermore  assume  that  each  component  of  the  spectrum  represents  an 
independent  random  process.  Fortunately,  for  the  physical  situation  we  are  addressing, 
each  component  does  represent  an  independent  source  of  turbulence;  5^^  (and  hence 
5^“^^)  represents  a  natural  cascade  of  turbulent  energy  from  long  to  short  wavelengths, 
whereas  (and  5'^*’^)  describe  the  turbulent  energy  distribution  among  modes  being 
driven  by  an  instability. 

3.2  Analysis  of  a  6-function  “Bump”  Component 

With  the  assumption  of  the  statistical  independence  of  the  “bump”  component 
of  the  spectrum,  we  now  proceed  to  analyze  its  effect  on  laser  propagation.  We 
begin  with  an  example  of  the  simplest  type  of  “bump”  in  the  turbulence  spectrum: 
a  6-function  “spike”  at  a  single  wavenumber  kb 

Sf^^(k}s2wSN^^S(k-kb)  (37) 

Here,  as  we  again  assume  two-dimensional  turbulence,  k  =  yjk^  +  k^,  and  the  factor 
of  2Tr/k  is  appropriate  for  normalization.  This  can  be  seen  by  inserting  (37)  into 
(8)  to  obtain  the  corresponding  spatial  correlation  function  of  the  refractive  index 
fluctuations 

R^ffis)  =  SN^Mkbs)  (.38) 
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with  the  level  of  the  fluctuations  given  by  (SN^)  =  SNq  (using  (14)). 

In  order  to  construct  the  phase  spectrum  induced  by  this  turbulent  spectrum  we 
first  use  the  exact  relation  (12)  to  find 


^2^2  y  ^  27r  SN^  i  6{k  -  h) 


/  sin 

V  k,zl2  ) 


The  integral  is  trivially  computed  to  give 


2klz^6N^  /sin  ^kl-klzl2\ 
-  kl  [  ^kl  -  klzl2  ) 


m 


(40) 


A  plot  of  is  shown  in  Fig. 3,  where  the  parameters  used  are  ko  =  6.28cm“*  (car¬ 
rier  wavelength  of  1cm),  2  =  10m  (propagation  distance  of  10m),  and  kb  =  6.28m“‘ 
(fluctuations  at  Im  wavelength).  Although  these  parameter  values  are  not  of  inter¬ 
est  for  laser  propagation  over  long  distances,  the  plot  serves  to  illustrate  the  general 
features  of  this  phase  spectrum.  It  is  obvious  that,  as  one  may  have  expected,  the 
phase  fluctuations  are  greatest  for  wavenumbers  near,  but  below  kb,  the  wavenum¬ 
ber  of  the  turbulent  fluctuations.  For  kx  ~  kb,  the  phase  spectrum  is  singular  (as  is 
S]^\k  =  kb)),  while  for  kx  >  kb  there  are  no  phase  fluctuations  at  all. 

Although  5^**^  is  non-negative  for  all  kx,  it  oscillates  with  increasing  frequency  as 
kx  approaches  kb  from  below.  The  zeroes  kn  of  Sf\kx)  are  easily  determined  to  be 

=  o<"<^  (41) 


which  have  a  separation  of  approximately  AA;„  ss  j {kbZ^).  Thus,  as  kbZ  becomes 

large,  the  oscillations  in  become  more  rapid.  This  has  the  following  consequence: 
given  kb  and  z,  the  value  of  at  0  <  fej,  <  fc;,  can  be  computed  from  (40).  However, 
for  large  kbZ  (as  is  our  usual  case  of  interest),  the  value  of  at  a  fixed  kx  is  very 
sensitive  to  the  value  of  2.  Indeed,  changing  2  by  an  amount  of  roughly  27r/A:i  can 
change  the  value  of  5^*”'  from  a  maximum  of  Sk^SN^Kk^  —  down  to  zero.  In 

practice,  this  means  that  changing  a  propagation  distance  of  tens  to  thousands  of 
kilometers  by  as  little  as  a  meter  could  change  the  phase  fluctuation  spectrum  at  any 
given  wavenumber  kx  <  kb. 

The  presence  of  these  rapid  oscillations  in  Sf\  together  with  the  fact  that  the 
phase  fluctuations  become  singular  at  kx  =  h  suggest  that  the  only  important  fluctua¬ 
tions  are  those  at  kx  =  kb.  Thus,  it  appears  we  are  justified  in  using  the  approximation 
(13)  for  the  phase  spectrum  induced  by  the  ^-function  “bump”  or 

Si'\kx)  =  klzSN{kx,kx  =  0)  =  2Kklz6N^j^  6{\kx\  -  kb)  (42) 

\f^x\ 
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Figure  3.  Plot  of  the  exact  phase  spectrum  (40)  induced  by  a  6-function  in 
the  turbulence  spectrum. 
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Even  though  a  singular  phase  spectrum  (such  as  either  (40)  or  (42))  introduces  some 
degree  of  difficulty  into  the  numerical  phase-screen  method  (where  the  value  of 
is  to  be  used  as  a  variance  for  a  probability  distribution),  the  approximate  expression 
(42)  does  permit  theoretical  analysis.  Furthermore,  we  shall  see  in  the  next  Section 
that  an  analysis  of  the  effect  of  such  a  “spike”,  bcised  on  the  use  of  (42),  is  useful  in 
predicting  the  effect  of  a  narrow  but  finite  “bump”  in  the  turbulence  spectrum. 

Taking  (42)  for  the  phase  spectrum,  we  can  trivially  compute  the  phase  correlation 
function  to  be 


=  /  ^  ^Trfc^z  8Nl  S{\k^\  -  kb) 


"^^0^  C\t2 


kb 


6Nq  cos  kbSx 


Thus,  the  level  of  the  phase  fluctuations  induced  by  the  “bump”  is 

{<s?)  =  ^  SK 

and  by  (16),  the  integral  scalelength  of  the  turbulence  is 

/ 

■  IT 


(43) 


(44) 


(45) 


where  £(,  =  2Tr/ kb  is  the  size  of  the  turbulent  eddies  corresponding  to  the  “bump” 
wavenumber  kb.  As  one  would  expect,  the  infinitesimally  narrow  phase  spectrum 

(42)  corresponds  to  a  phase  correlation  function  which  has  infinite  range.  A  realistic 
“bump”  phase  spectrum  will  always  have  some  finite  width  AA:  and  will  therefore 
correspond  to  a  correlation  function  which  decays  over  a  distance  proportional  to 
{Ak)~^-,  as  long  as  the  propagation  distance  z  of  interest  is  much  greater  than  this 
decay  length,  we  are  justified  in  using  the  results  of  this  analysis  bcised  on  (42). 

To  investigate  the  effects  of  such  a  “bump”  on  the  propagation  of  a  laser  beam,  we 
consider  the  predictions  for  the  ensemble- averaged  wave  intensity  (/(i,z))  and  wave 
spectrum  {I{k,z))  given  by  (18)  and  (30)  of  Section  2.3,  respectively.  We  assume 
that  the  beam  has  an  initial  Gaussian  wave  profile  (as  given  by  (22))  so  that  the 
amplitude  correlation  Ma{s',z)  is  given  by  (23).  For  the  phase  correlation  function 

(43)  then,  the  ensemble- averaged  wave  intensity  is 
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Clearly,  this  result  reduces  to  the  vacuum  propagation  result  (25)  in  the  absence 
of  turbulence  (since  the  Bessel  functions  satisfy  /o(0)  =  1  and  /„(0)  =  0  for  n  ^  0). 
If  only  the  leading  term  (proportional  to  Iq)  were  present,  the  mean  influence  of 
the  turbulence  would  be  just  a  pure  damping  (by  a  factor  of  Ioi{<f>l))  exp{  — {4>l)))  of 
the  vacuum  propagation  result  (although  this  would  not  conserve  the  beam  power). 
Besides  damping  the  wave  intensity  near  the  beam  axis  (at  1  =  0),  however,  the  tur¬ 
bulence  causes  the  beam  to  spread  over  and  above  that  caused  by  natural  diffraction, 
and  it  develops  intensity  maxima  off-axis  at  the  transverse  locations  ±x„  given  by 

x„(^)  ^  f  (47) 

A  plot  of  this  result  is  shown  in  Fig.  4a,  where  it  is  compared  with  the  vacuum 
propagation  profile.  For  the  purpose  of  illustration,  we  have  chosen  the  following 
parameters 

ko  =  6280cm“^  (A  =  10//m) 

a  =  1.0m  =>  To  =  0.71m 

k\,  —  6.28m“®  (4  =  Im) 

2  =  1000km 

{<^1)  =  3.14  =>  =  2.5  X  10"*^  from  (44)  (48) 

With  these  values,  the  intensity  peaks  (47)  are  seen  to  be  spaced  at  five  meter  inter¬ 
vals.  With  the  use  of  Bessel  function  summation  identities,  one  can  easily  show  that 
the  integral  of  the  intensity  (46)  over  x  is  indeed  equal  to  its  value  at  z  =  0. 

The  prediction  for  the  ensemble- averaged  wave  spectrum  {Ii,{k,z))  produced  by 
the  turbulence  (30)  requires  the  Fourier  transform  (31)  of  the  turbulence  modulation 
Mt  or 

Mrik-z)  = 
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Figure  4.  Plot  of  the  predicted  ensemble- averaged  intensity  produced  by 
a  (5-function  in  the  turbulence  spectrum,  with  parameters  given  by  (48). 
Vacuum  propagation  result  shown  by  lighter  curve,  (a)  Spatial  intensity, 
(b)  Wave  spectrum. 
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/oo 

ds  /o((<^^))  +  2  X:  In{{d>l))cos{nhs) 

n=l 

=  2ie-W>{/„({^J))S(t) 

OO 

+  E  WD)  +  Sik  +  «A:0]}  m 


n=l 


With  the  Fourier  transform  ip{k,0)  of  the  initial  Gaussian  profile  (22) 


(50) 


the  form  of  (49)  permits  the  trivial  evaluation  of  (30)  to  give 
{hik,z)}  =  27ra^e-<'>'l>{loi{d>l})e-'‘"<^" 

OO 

+  E  WD)  (51) 


n=l 


Again,  this  result  reduces  to  the  vacuum  expression  (the  square  of  (50))  in  the  absence 
of  turbulence,  and  the  leading  term  damps  this  value  in  the  presence  of  turbulence. 
Like  the  spatial  intensity  (46),  the  wave  spectrum  is  also  spread  by  the  turbulence, 
and  develops  sidebands  at  harmonics  dckn  of  the  turbulence  “bump”  wavenumber  ki,: 


(52) 


A  plot  of  this  result  is  shown  in  Fig.  4b,  for  the  same  parameters  (48). 

Finally,  one  can  compute  from  (46)  the  evolution  of  the  width  of  the  beam  as  it 
propagates  As  defined  in  (26)  we  have 

>■^(2)  =  J 

=  ^+2e-<«£x;/„((«)  (53) 

n=l 

Here  we  have  computed  the  Gaussian  integrals  and  have  again  used  the  Bessel  function 
summation  rule.  Since  by  (47)  the  intensity  peaks  at  ±Xn  are  proportional  to  n,  we 
use  the  Bessel  function  generating  function  to  find 


^  n^/„(<)  =  |te‘ 


so  that  we  have 


r^(z) 


n=:l 


2  V  /  4/fc2  2 


(T‘ 

y 


1  + 


k^a^ 


+ 


(54) 


(55) 


where  we  have  used  the  relation  (44)  between  and  Identifying  from  (55)  the 
initial  width  tq  =  cr/\/2,  this  expression  is  exactly  that  given  in  (27)  (indicating  the 
cubic  growth  of  the  beam  width  due  to  turbulence)  with  the  turbulent  broadening 
length  Lt  given  by 


Lt  = 


■  ■ 

1/3 

kbSNl 

tSN^ 

1/3 


(56) 


VVe  note  that  this  definition  coincides  precisely  with  the  definition  (28);  it  is  trivial 
to  compute  from  (37)  that  the  transverse  scalelength  it  (defined  by  (29))  is  exactly 
4,  and  from  (45)  we  have  £,  =  ^/tt.  For  the  parameters  (48)  above,  Lj  =  234km 
while  the  Rayleigh  length  is  Lr  =  kor^  =  314km;  thus  in  this  case,  it  is  clear  that  the 
turbulent  broadening  dominates  over  natural  diffraction,  as  borne  out  by  Fig.  4. 


3.3  Numerical  Simulation  for  a  Gaussian  “Bump” 

A  numerical  simulation  of  propagation  through  a  turbulent  medium  characterized 
by  the  singular  spectral  density  (37)  (or  corresponding  singular  phase  spectrum  (40) 
or  (42))  using  the  phase-screen  technique  is  difficult  in  that  it  utilizes  fast  Fourier 
transform  methods.  Therefore,  we  have  considered  a  spectral  density  “bump”  of  the 
form 

5g)(fc)  =  (57) 

KfCl) 

with  the  hope  that  if  the  width  k  is  small  compared  with  the  “bump”  wavenumber  kf, 
then  the  results  of  the  previous  Section  will  apply.  We  therefore  assume  the  following 
to  hold 


kb  =  2TTjib  K  =  27r/Z/o 

h  ^  K  =>  ib  (58) 

where  the  inverse  of  the  A:-space  width  of  the  spectral  density  gives  the  range  (or 
outer  scale)  Lq  of  the  spatial  correlation  function.  In  the  numerical  simulations,  we 
will  use  the  values 

ib  =  Im  (A:6  =  6.28m“*) 

Lo  =  30m  («:  =  0.2lm~^)  (59) 

A  plot  of  5^^  with  these  parameters  is  shown  in  Fig.  5a. 

In  the  regime  (58),  the  two-dimensional  Fourier  transform  (8)  can  be  approximated 
analytically  to  give 

R^j^Hs)^SN^Jo{kbs)  (60) 

which  is  the  same  as  that  for  the  6-function  “bump”  (38).  The  exact  transform 
(computed  ninnerically)  shown  in  Fig.  5c  behaves  like  (60)  for  small  s,  but  is  observed 
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Figure  5.  Turbulence  functions  used  in  numerical  simulations,  (a) 
from  (57)  with  parameters  (59).  (b)  Log-log  plot  of  (c)  Spatial 

correlation  function  /?j^^(s).  (d)  Phase  correlation  function  R^^\sx)- 


to  decay  mu<'h  faster  than  Jo\  indeed,  the  correlation  function  in  this  case  should 
exhibit  a  much  shorter  range  (~  Lq)  than  the  correlation  function  corresponding  to 
a  (5-function,  due  to  the  finit(‘  width  of  the  spectral  density. 

As  we  have  justified  in  the  preceding  Sections,  we  take  the  spectral  density  for  the 
phase  fluctuations  to  be  approximated  by  (13)  or 

sf{k,)  =  klzSN{K,K  =  Q)  =  (61) 

for  which  the  phase  correlation  function  is 


R^^\s,)  =  SN^  cos  hsr  (62) 

fCf) 

This  is  plotted  in  Fig.  5d.  Again,  this  form  is  similar  to  f43)  for  the  (5-function  spectral 
density,  but  now  explicitly  exhibits  a  correlation  decay  over  the  length  Lq  =  27r//c. 
The  level  of  phase  fluctuations  ((j5j)  and  the  integral  scalelength  ii  are  the  same  as 
those  (44,45)  fo'-  the  ^-function 

(fi|)  =  <>(0)  =  ‘^6Ni  ^  «,  =  !  =  -  (63) 

while  the  transverse  scalelength  it  (from  (29)  using  (57))  is 

l,  =  » £,  (64) 


for  To  >  4- 

As  in  the  previous  Section,  we  would  like  to  be  able  to  construct  the  predic¬ 
tions  (46)  and  (51)  for  the  ensemble-averaged  wave  intensity  and  wave  spectrum  for 
this  type  of  spectral  “bump”.  Even  though  this  is  not  possible  analytically,  we  can 
construct  these  quantities  numerically.  Furthermore,  with  the  non-singular  spectral 
density  (61)  we  can  use  the  phase-screen  method  to  numerically  simulate  the  prop¬ 
agation  of  a  laser  beam  through  a  realization  of  the  turbulence  described  by  (57). 
Simulating  the  propagation  of  a  number  of  shots  through  an  ensemble  of  realizations 
of  the  turbulence  and  averaging  over  the  final  intensities  and  wave  spectra  will  al¬ 
low  us  to  check  the  predictions  based  on  a  numerical  construction  of  (7(,(.'c,r))  and 
{h{k,z))  from  (62),  as  well  as  the  predictions  given  in  the  previous  Section  based  on 
a  ^-function  spectral  “bump”. 

Our  numerical  results  are  shown  in  Fig.  6.  Using  the  parameters  given  in  (48,59), 
the  intensity  of  a  beam  propagated  through  a  single  realization  of  the  turbulence  is 
shown  in  Fig.  6a.  The  wave  intensity  clearly  shows  the  effect  of  the  loss  of  coher¬ 
ence  due  the  turbulence;  it  has  developed  “spikes”  (or  “beamlets”)  and  has  spread 
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Figure  6.  Numerical  simulation  of  propagation  through  a  single  realization 
of  Gaussian  “bump”  turbulence,  with  parameters  given  by  (48,59),  compared 
with  prediction  based  on  (62)  (smooth  dark  curve)  and  vacuum  result  (light 
curve),  (a)  Intensity,  (b)  Wave  spectrum. 
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much  more  than  the  vacuum  result  (shown  by  the  lighter  curve).  'Flie  prediction  for 
the  ensemble- averaged  intensity  (shown  by  the  smooth  dark  curve)  was  numerically 
constructed  from  (18)  using  the  Gaussian  “bump”  phase  correlation  function  (62); 
this  prediction  compares  quite  favorably  with  that  shown  in  Fig.  4a  based  on  the 
5-function  “bump”.  The  wave  spectrum  after  propagation  through  a  single  turbu¬ 
lence  realization  (shown  if  Fig.  6b)  has  clearly  developed  maxima  at  harmonics  of 
the  “bump”  wavenumber  fcj,  and  again,  the  ensemble- average  prediction  numerically 
constructed  for  the  Gaussian  “bump”  case  is  very  similar  to  the  analytical  form  for 
the  5-function  case  shown  in  Fig.  4b. 

An  ensemble  average  over  100  shots  (propagation  through  100  realizations  of  the 
turbulence,  averaged  numerically)  is  shown  in  Fig.  7.  Now  the  numerical  ensemble 
average  is  much  smoother  than  the  single-reahzation  result  of  Fig.  6  (the  “spikes”  or 
“beamlets”  have  been  smoothed  out  by  averaging,  as  to  be  expected),  but  the  effect 
of  the  turbulence  producing  increased  beam  spreading  is  still  evident.  However,  it 
is  also  interesting  that  the  off-axis  intensity  maxima,  as  predicted  by  both  the  5- 
function  and  Gaussian  “bump”  model  for  the  ensemble- averaged  intensity,  do  not 
appear  in  the  actual  numerical  results.  The  numerical  ensemble- average  does  spread 
just  as  predicted,  but  the  result  appears  to  wash  out  the  off-axis  maxima  with  a  sort 
of  average  behavior.  The  wave  spectrum,  on  the  other  hand,  very  closely  matches  the 
prediction  for  maxima  at  harmonics  of  the  “bump”  wavenumber  k^. 

The  discrepency  between  the  numerical  and  predicted  results  can  be  explained  as 
follows:  The  prediction  (18)  is  based  on  a  theoretical  analysis  of  a  single  phase-screen 
step.  That  is,  one  can  write  down  a  closed-form  expression  for  the  final  intensity  (say, 
at  2  =  1000km  as  in  our  example)  as  if  it  were  propagated  in  vacuum  for  the  first  half 
of  the  distance  (t.e.,  500km),  multiplied  by  a  random  phase  function  exp(z^(x;2)) 
which  accounts  for  the  turbulence  over  the  entire  distance  of  propagation  z,  and  then 
propagated  in  vacuum  the  remaining  distance  (another  500km).  Another  prediction 
could  be  generated  by  first  multiplying  by  the  random  phase  function  (at  2  =  0)  and 
then  propagating  in  vacuum  the  entire  distance  (2  =  1000km).  This  latter  analy¬ 
sis  would  have  produced  an  expression  for  the  predicted  ensemble- averaged  intensity 
{l{x,z))  similar  to  that  in  (18),  except  for  the  factor  of  1/2  in  the  argument  of  Mj 
would  be  unity.  Using  this  alternative  prediction  in  the  5-function  “bump”  analy¬ 
sis  of  Section  3.2  would  yield  off-axis  intensity  maxima  at  =  nkbzfko  (t.e.,  at 
±10,  ±20  . . .  meters  for  our  parameters),  instead  of  at  =  nkb{zf2)lkQ  (±5,  ±10  . . . 
meters)  as  in  (47).  Of  course,  in  order  to  conserve  beam  power,  the  intensity  at  each 
of  these  maxima  would  be  lower  since  the  profile  is  more  spread  out.  With  this  “ini¬ 
tial”  phase-screen  approach,  however,  the  expression  (30)  for  the  wave  spectrum  is 
unchanged;  thus,  the  wave  spectrum  produced  by  “bump”  turbulence  develops  max¬ 
ima  at  harmonics  of  the  “bump”  wavenumber  in  both  analyses.  The  interpretation 
is  quite  clear  then:  a  single  application  of  the  phase-screen  representing  a  spectral 
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Figure  7.  Numerical  ensemble  average  of  100  shots  through  different  real¬ 
izations  of  Gaussian  “bump”  turbulence,  with  parameters  given  by  (48,59), 
compared  with  prediction  based  on  (62)  (smooth  dark  curve)  and  vacuum 
result  (light  curve),  (a)  Intensity,  (b)  Wave  spectrum. 
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bump  at  ^  =  kh  (irregardless  of  where  it  is  placed)  will  cause  the  wave  spectrum  to 
develop  peaks  at  harmonics  ±A;„  =  ±.nki,.  This  can  be  thought  of  as  a  scattering 
of  rays  to  a  transverse  wavenumber  of  A:„;  these  rays  then  propagate  with  an  angle 
=  kn/ko  —  nkbiko  to  the  2-axis.  Now,  if  the  phase-screen  is  applied  at  the  midpoint 
(z/2),  then  the  transverse  distance  reached  by  these  rays  is  Xn  =  0nzl2  =  nkf,zf{2kQ) 
(producing  intensity  maxima  as  predicted  by  (47)).  If  the  phase-screen  is  applied  at 
2  =  0,  however,  the  rays  propagate  twice  as  far  at  the  new  angles  this  results  in 
intensity  maxima  at  the  more  widely  spaced  locations  of  Xn  =  nkbz/ko- 

Naturally,  in  reality,  propagation  through  turbulence  is  a  continuous  process,  not 
accurately  modeled  by  the  imposition  of  a  single  random  phase- shifting  screen.  Nu¬ 
merically,  we  simulate  the  propagation  by  using  a  series  of  intermediate  phase-screens, 
with  the  idea  that  as  the  number  of  phase-screens  increases  (with  decreasing  separa¬ 
tion),  the  more  closely  the  effect  of  continuous  propagation  is  simulated.  Note  that 
since  by  (16)  the  level  of  the  phase  fluctuations  is  proportional  to  the  propagation 
distance  2,  it  is  possible  to  divide  the  total  distance  2  up  into  N  phase-screen  inter¬ 
vals  of  length  z/N  (with  a  level  of  the  phase  fluctuations  on  each  screen  proportional 
to  z/N)  so  that  the  total  level  of  the  phase  fluctuations  remains  constant  as  N  is 
increased.  In  this  way,  the  influence  of  the  turbulence  on  the  beam  (or  the  scattering 
of  the  rays)  becomes  more  gradual.  The  goal  in  the  numerical  phase-screen  method 
then  is  to  find  an  optimal  separation  distance  for  the  phase-screens;  i.e.,  we  want  to 
minimize  the  number  of  phase-screens  required  (to  decrease  the  computation  time) 
with  the  property  that  increasing  the  number  does  not  significantly  change  the  result. 
In  fact,  in  the  parameter  regime  we  are  exploring,  we  have  found  that  phase-screens 
separated  by  about  100km  is  sufficient,  with  (according  to  (48))  a  phase  fluctua¬ 
tion  level  of  about  «  0.3  on  each  phase-screen;  that  is  the  number  used  in  the 
simulation  results  reported  in  this  article. 

So  now  the  discrepancy  between  the  simulation  results  and  the  predicted  ensemble- 
averaged  intensity  is  clear:  at  each  of  the  ten  phase-screens  (located  at  2^),  rays 
are  scattered  to  angles  d„  =  nkb/ko,  thence  to  propagate  the  remaining  distance 
2  —  2,  and  produce  an  intensity  peak  at  (2  —  z,)^„  =  (2  —  Zi)nkb/ko.  Since  there  are 
ten  values  of  Zi,  there  are  ten  intensity  peaks  associated  with  each  value  of  n:  this 
effectively  smears  out  the  off-axis  intensity  maxima.  If  the  number  of  phase-screens 
were  increased  (be.,  more  closely  approximating  the  actual  physical  propagation), 
the  initensity  peaks  would  become  more  closely  spaced,  effectively  becoming  smeared 
so  that  indeed  no  peaks  are  present.  Thus,  it  appears  that  using  ten  phase  screens 
is  enough  to  simulate  the  true  continuous  propagation  (the  smearing  appears  to  be 
complete).  Furthermore,  we  note  that  while  the  prediction  for  the  ensemble- averaged 
intensity  (based  on  the  single  mid-point  phase-screen  analysis)  appears  to  be  quite 
accurate  (in  a  mean  sense,  averaging  out  the  off-axis  intensity  maxima  which  are  not 
physical);  indeed,  our  research  has  shown  that  this  prediction  is  much  better  than 
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one  based  on  a  single  phase-screen  at  ^  =  0.  Perhaps  more  importantly,  even  though 
the  ofT-axis  intensity  maxima  do  not  occur,  the  prediction  for  the  growth  of  the  beam 
width  in  terms  of  Lr  (depending  on  the  “bump”  parameters  £,  =  ^/tt  and  it  ^  ib) 
appears  to  be  borne  out  by  the  simulations. 

3.4  A  “Bump”  Superimposed  on  a  Cascade 

In  th  is  Section  we  will  now  consider  the  case  of  a  turbulence  spectral  density 
composed  of  both  a  cascade  and  a  “bump”  component.  We  therefore  begin  with  a 
spectral  density  Sjw  of  the  form 

SN{k)  =  Sffik)-{-S]^\k)  (65) 

for  which  the  level  of  refractive  index  fluctuations  (or  total  power)  is 

r  fPk 

Since  the  effects  on  beam  propagation  (such  as  turbulent  broadening)  depend  on  the 
magnitude  of  {SN'^),  it  is  clear  from  (66)  that  increasing  the  power  in  the  “bump” 
component  of  the  spectrum  (i.e.,  increasing  {SN^))  while  holding  the  cascade  compo¬ 
nent  fixed  will  cause  the  beam  to  be  more  severely  affected  by  the  turbulence.  For 
example,  increasing  {SN^)  (and  thus  the  total  {SN'^))  will,  by  (16),  increase  the  level 
of  phase  fluctuations  {(f>‘^)  and,  by  (28),  decrease  the  turbulent  broadening  length  Lt 
so  that  the  spreading  caused  by  the  turbulence  will  begin  to  occur  over  a  shorter 
distance.  Furthermore,  as  is  increeised  to  levels  much  higher  than  {SN^),  it  is 

reasonable  to  expect  that  the  effect  on  the  beam  will  be  dominated  by  the  turbulence 
in  the  “bump”  component.  Thus,  for  {SN^)  ^  ^^e  integral  and  transverse 

scalelengths  (which  also  affect  both  {(f>^)  and  Lt)  will  become  more  representative 
of  the  “bump”  component;  for  a  Gaussian  “bump”  as  in  the  preceding  Section,  we 
would  expect  ii—^ibl'K  and  it—^ib-  In  this  case  then,  for  a  spectral  “bump”  at  large  k 
(or  short  wavelength)  as  depicted  in  Fig.  1,  we  would  expect  these  scalelengths  to  be 
smaller  than  if  the  cascade  component  were  dominant;  thus,  Lj  is  further  decreased 
and  turbulent  broadening  becomes  a  more  severe  effect.  Even  though  we  have  seen 
in  the  previous  Section  that  off-axis  intensity  maxima  (as  predicted  using  the  single 
midpoint  pha.se-screen  theory)  do  not  occur,  the  prediction  for  the  cubic  growth  of 
the  beam  width  in  terms  of  Lj  (depending  on  the  “bump”  parameters  it  and  it)  still 
holds. 

In  order  to  directly  compare  the  relative  effects  of  the  cascade  component  5^* 
and  the  “bump”  component  we  can  vary  the  magnitude  of  each  component 

w'hilc  holding  the  total  power  in  the  turbulence  spectrum  (or  {SN^))  fixed.  We  thus 
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(67) 


define  the  fraction  q  of  the  total  power  in  the  “hump”  component  to  he 

'  “  (Sm)  (SN^) 

SO  that  {8N'^)  by  (66)  is  constant  as  q  is  varied.  Now  since  (65)  is  a  linear  superposition 
of  cascade  and  “bump”  components,  so  is  the  phase  spectrum  produced  by  (65). 
Therefore,  it  follov/s  that  the  level  of  phase  fluctuations  is 

{4?)  =  klz  SN(k)  =  kiz{SN^)t, 

=  +  ime;] 

=  +  (1  - 

=  kU{SN^){t‘  -  ,(f?  -  «*)]  (68) 


As  mentioned  above,  if  the  “bump”  component  is  at  larger  wavenumber  than  the 
cascade  (as  depicted  in  Fig.  1),  then  the  scalelengths  characterizing  the  “bump” 
component  (such  as  and  i^^)  will  be  smaller  than  those  of  the  cascade  component. 
Thus,  we  see  from  (68)  that  as  power  is  transferred  into  the  bump  component  (q  is 
increased),  the  phase  fluctuations  measured  by  ((^^)  are  decreased  (for  >  i^).  This 
was  previously  noted  by  Goldring. ^ 

While  this  result  seems  to  indicate  that  taking  energy  out  of  the  cascade  and 
putting  it  into  the  “bump”  tends  to  diminish  the  overall  effect  of  the  turbulence,  one 
must  remember  that  even  though  this  transfer  of  energy  decreaises  the  magnitude  of 
the  phase  fluctuations,  the  fluctuations  are  of  shorter  wavelength.  As  the  scalelength 
of  the  phase  fluctuations  becomes  smaller  (especially,  small<=*r  than  the  beam  width), 
the  loss  of  transverse  coherence  of  the  beam  becomes  more  rapid.  The  combination  of 
these  effects  is  embodied  in  the  definition  (28)  of  the  turbulent  broadening  length,  Lt, 
for  which  one  also  needs  to  compute  the  transverse  scalelength  it  of  the  turbulence. 
Thus,  while  the  effective  integral  scalelength  it  is  given  by  (68)  to  be 

?,  =?<?  +  (! (39) 


the  effective  transverse  scalelength  is  from  (29) 

V9"/  jdkk^Srj(k)  J dk  k^  {Sli\k)  +  sii\k)\ 

Using  the  first  relation  to  define  the  transverse  scalelength  of  each  component  .sepa¬ 
rately  this  becomes 


i, 


{8N^) 
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1+9 


from  (28)  the  turbulent  broadening  length  is 


7r2(^Ar2) 


(^rr  1+9 


(£nv^^ 


where  and  Lj  are  defined  in  terms  of  the  total  power  {SN'^).  Therefore  we  see 
that  if  the  characteristic  scalelength  of  the  cascade  is  greater  than  that  for  the 
“bump”,  then  increasing  the  power  in  the  “bump”  (increasing  q)  will  decrease  the 
broadening  length  Lj- 

We  have  explored  this  effect  numerically  using  the  phase-screen  technique  to  prop¬ 
agate  a  beam  through  turbulence  described  by  a  composite  spectral  density  such  as 
(65).  The  “bump”  component  5^^  was  taken  to  be  the  Gaussian  “bump”  discussed 
in  the  previous  Section,  while  the  cascade  component  was  modeled  by  a  von  Karman 
spectrum 

^  [1  +  (k/kminYY/^ 

In  the  regime  kmin  <C  A;  <C  kmax,  this  spectrum  behaves  like  the  two-dimensional 
Kolmogorov  cascade  5^*(fc)  ~  The  inner  scale  /§  and  outer  scale  Lq  of  the 

cascade  were  taken  to  be 

/q  —  27r  j  kfTfiax  —  0.5m 

LI  =  2'Klkmin  =  10m  (74) 

The  cascade  spectrum  and  its  corresponding  correlation  functions  are  shown  in  Fig.  8. 
The  turbulence  scalelengths  for  this  cascade  were  computed  numerically  to  be 


—  5.1m 

=  9.0m  (75) 

For  the  Gaussian  “bump”  parameters  (59),  the  turbulence  scalelengths  from  (63,64) 


£*  =  0.3m 
=  1.0m 

(4)7#;  =  3.3m 
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Figure  8.  Cascade  model  used  in  numerical  simulations,  (a)  S\^\k)  from 
(73)  with  parameters  (74).  (b)  Lug-log  plot  of  S\^\k).  (c)  Spatial  correlation 
function  R]^\s).  (d)  Pha.se  correlation  function 
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With  these  values,  the  turbulent  broadening  length  for  the  pure  cascade  spectrum 
(q  =  0)  is  Lj  =  330km,  while  for  the  pure  “bump”  spectrum  {q  =  1)  it  is  =  231km. 
Thus  in  this  case,  we  should  ob.serve  a  more  pronounced  broadening  of  the  beam  as 
energy  is  transferred  from  the  cascade  to  the  “bump”. 

The  beam  with  the  parameters  given  in  (48)  was  propagated  using  the  phase- 
screen  method  as  described  in  previous  Sections,  with  phase-screens  located  every 
100km  for  a  1000km  total  propagation  distance.  We  present  here  the  results  of  nu¬ 
merically  averaging  over  an  ensemble  of  100  shots  through  different  realizations  of  the 
turbulence  at  fixed  q.  The  numerical  ensemble-averaged  intensity  and  wave  spectrum 
are  compared  with  predictions  constructed  numerically,  by  (18)  and  (30),  using  the 
actual  phase  correlation  function  for  each  q. 

The  effect  of  just  the  pure  cascade  spectrum  (q  =  0)  is  shown  in  Fig.  9.  Obviously, 
the  beam  has  been  broadened  (in  both  x-space  and  fc-space)  due  to  the  turbulence, 
but  not  to  the  degree  that  it  was  broadened  by  the  pure  “bump”  turbulence  (^  =  1) 
as  shown  in  Fig.  7.  This  is  consistent  with  Lt  for  the  pure  cascade  being  larger 
than  that  for  the  pure  “bump”  spectrum.  Furthermore,  the  wave  spectrum  in  this 
case  does  not  exhibit  the  harmonic  structure  characterizing  that  for  the  “bump” 
case;  this  is  presumably  because  the  cascade  spectrum  is  centered  at  A:  =  0  and  is 
much  broader  than  the  “bump”  spectrum.  Finally,  it  is  apparent  that  the  expressions 
for  the  ensemble- averaged  intensity  (18)  and  wave  spectrum  (30),  based  on  a  single 
midpoint  phase-screen  theory,  give  very  accurate  predictions  for  these  quantities. 

In  Fig.  10  we  show  the  composite  spectral  densities  for  q  =  0.325,  0.7,  and  0.83. 
These  values  were  chosen  because  the  magnitude  of  the  “bump”  component  relative 
to  that  of  the  cascade  component  (be.,  S^^\kb) / for  these  cases  is  0.1,  0.5  and 
1.0,  respectively.  The  results  of  propagating  through  each  of  these  turbulence  models 
are  shown  in  Figs.  11-13.  The  trend  toward  a  broader  beam  as  q  increases  is  clear, 
as  predicted  above.  At  (?  =  0.7,  the  wave  spectrum  begins  to  show  the  harmonic 
structure  of  the  pure  “bump”  case;  the  prediction  for  the  spatial  intensity  also  begins 
to  develop  off-axis  maxima,  although  the  numerical  propagation  results  show  that 
these  should  not  be  present  (as  explained  in  the  preceding  Section).  By  q  =  0.83,  the 
beam  is  broadened  almost  as  much  as  in  the  pure  “bump”  case;  indeed,  by  (72),  the 
value  of  Lt  is  245km  for  this  case,  compared  with  234km  for  the  pure  “bump" . 
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Figure  9.  Numerical  ensemble  average  of  100  shots  through  different  real¬ 
izations  of  cascade  turbulence  (73,74)  and  Fig.  8,  compared  with  prediction 
based  on  (18,30)  (smooth  dark  curve)  and  vacuum  result  (light  curve). 

(a)  Intensity,  (b)  Wave  spectrum. 
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Figure  10.  ( 
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Figure  11.  Numerical  ensemble  average  of  100  shots  through  different  re¬ 
alizations  of  turbulence  with  composite  spectral  density  for  q  =  0.325. 

(a)  Intensity,  (b)  Wave  spectrum. 
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Figure  12.  Numerical  ensemble  average  of  100  shots  through  different  re¬ 
alizations  of  turbulence  with  composite  spectral  density  for  q  =  0.7. 

(a)  Intensity,  (b)  Wave  spectrum. 
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Figure  13.  Numerical  ensemble  average  of  100  shots  through  different  re¬ 
alizations  of  turbulence  with  composite  spectral  density  for  q  =  0.83. 

(a)  Intensity,  (b)  Wave  spectrum. 
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3.5  The  Effect  of  a  “Bump”  in  Three  Dimensions 

Throughout  this  report  we  have  focused  for  simplicity  on  the  effects  of  '‘bump" 
turbulence  on  propagation  in  two  dimensions.  Based  on  the  analysis  in  the  yjievious 
Section,  we  will  now  comment  on  the  effect  of  a  “hump”  superimposed  on  a  cascade 
spectrum  in  three  dimensions.  For  the  purpo.se  of  analytical  computation,  we  take 
the  isotropic  cascade  spectrum  to  be  a  Gaussian 

=  (77) 

while  the  “bump”  is  also  a  Gaussian  centered  at  the  “bump”  wavenumber  kt, 

=  (78) 

In  these  expressions,  A:  ~  0  is  the  radius  in  the  three-dimensional  k-space.  As  in  the 
previous  Section,  the  total  turbulence  spectrum  is  simply  the  sum  of  5^'  and  5^^ 
The  total  power  in  each  component  of  the  spectrum  can  be  computed  to  be 

(80) 

which  is  the  three-dimensional  analogy  of  (66).  Here  we  have  also  made  the  assump¬ 
tion  that  the  “bump”  component  is  very  narrow  and  located  at  large  wavenumber, 
kb  >>  Kb  (just  as  in  (58)).  At  this  point  it  is  interesting  to  compare  the  relative  power 
in  the  two  components; 

^  -  (SNJ)  SJ 

Consider  the  case  where  the  “bump”  component  represents  turbulence  at  a  much 
smaller  scalelength  than  the  cascade;  for  example,  let  us  take  the  following  parame¬ 
ters: 


kb  =  2-k  jib  =  27r/(lm)  =  6.28m“* 

Kb  =  2-k/ Lb  =  27r/(10m)  =  0.63m“^  (82) 

Kc  =  27r/Lc  =  27r/(10km)  =  6.28  X  10~‘’m“' 


Thus  the  turbulence  can  be  thought  of  as  turbulent  eddies  in  the  range  of  G  =  Im 
to  Lb  =  10m  (produced  by  some  instability,  for  example)  superimposed  on  a  cascade 
of  turbulence  with  scalelength  of  Lc  =  10km  (produced  by  some  other  large-scale 
mechanism).  In  this  case,  the  relative  power  (81)  becomes 


p  =  2  X  10" 


5? 


(S3) 
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'riuis,  even  if  the  magnitude  5°  =  S]^\ki,)  of  the  “bump”  component  is  a  factor  of 
10‘°  less  than  that  of  the  cascade  =  6'^^^(0),  the  “bump”  component  still  has  20 
times  the  power  of  the  cascade  component.  An  example  of  this  spectrum,  with  the 
parameters  as  in  (82),  is  shown  in  F'ig.  14.  In  this  case,  the  fraction  </  of  the  total 
power  [SN'^)  =  (SN^)  +  (SN^)  that  resides  in  the  “bump”  component  is 


<?  = 


P 

1  +  p 


=  0.95 


(841 


This  is  a  purely  geometric  effect  of  the  three-dimensional  k-space:  while  the  “bump” 
is  narrow  and  of  much  smaller  magnitude,  it  is  located  at  large  radius  in  k-space  and 
thus  the  volume  factor  works  to  increa.se  the  power  in  the  “bump”  component. 

The  effect  of  the  “bump”  component  can  now  be  assessed  in  terms  of  the  turbulent 
broadening  length  Tj,  using  (72).  For  this  we  need  the  characteristh’  scalelengths 

=  4.0  X  lO^m  4/2  =  0.5m 

4  =  Tc/\/2  =  7.0  X  lO^m  =  4  =  TOm  (85) 

=>  (4)74  =  =  1.25  X  lO^m  =J>  (4)74  =  24  =  2.0m 


which  are  computed  from  the  spectral  densities  using  the  three-dimensional  versions 
of  (16)  and  (29).  These  values  give 


(4)V4  _  1 
(4)74  2V2«6 


=  6.3  X  10=^ 


(86) 


Together  with  the  value  of  q  =  0.95  for  this  example,  we  see  that  the  ratio  of  the 
relevant  scalelengths  is  such  that  the  turbulent  broadening  length  (72)  is  decreased  by 
a  factor  of  almost  twenty,  from  the  value  characteristic  of  the  cascade  component 
to  the  value  determined  by  the  “bump”  component.  Thus,  in  this  case  where  the 
“bump”  component  appears  to  be  an  insignificant  contribution  to  the  total  spectrum 
(as  in  Fig.  14),  not  only  does  almost  all  of  the  power  reside  in  the  “bump”,  but  it 
completely  determines  the  effects  of  turbulence  on  the  beam;  the  “large”  cascade 
component  can  be  ignored. 
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Figure  14.  A  three-dimensional  isotropic  spectrum  {k  =  |k|)  with  a  “bump” 
superimposed  on  a  cascade  (given  by  (77)  and  (78)).  The  parameters  are 
those  in  (82),  with  5°  =  10^°m^  and  5°  =  Im^.  In  this  case,  95%  of  the 
power  in  the  spectrum  resides  in  the  “bump”  component. 
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4  Conclusions 

In  this  report  we  have  anlayzed  the  effect  of  a  “bump” -type  turbulence  spectrum 
on  laser  propagation,  and  we  have  compared  our  results  with  those  produced  by  the 
standard  assumption  of  a  cascade-type  spectrum.  In  general,  the  effect  of  narrow 
“bump”  in  the  turbulence  spectrum  (at  a  finite  wave  number  kb)  is  to  generate  side¬ 
bands  in  the  wave  spectrum  at  multiples  ±nkb  of  the  “bump”  wavenumber.  On  aver¬ 
age,  the  spatial  broadening  of  the  beam  behaves  as  predicted  (by  a  theory  based  on  a 
modified  Huygens  principle)  and  spreads  with  the  cube  of  the  propagation  distance; 
the  characteristic  broadening  distance  Lt  for  this  effect  decreases  as  the  scalelength 
of  the  turbulent  eddies  {b  =  decreases  and  as  the  level  of  the  refractive  index 

fluctuations  {SN^)  increa.ses.  On  this  basis  we  can  conclude  that  given  a  cascade  spec¬ 
trum  and  a  “bump”  spectrum  with  the  same  level  of  fluctuations  {SN'^),  the  spectrum 
with  the  smaller  characteristic  scalelength  will  have  the  greater  effect  on  the  beam 
(because  Lj  for  that  spectrum  will  be  smaller,  and  hence  turbulence  effects  will  be 
observed  over  shorter  propagation  distances).  Since  by  definition  a  cascade  spectrum 
is  dominated  by  a  peak  at  /;  =  0  and  usually  decreases  rapidly  as  some  power  law 
at  large  k  (representing  a  cascade  of  energy  from  longer  to  shorter  wavelengths),  it 
is  likely  that  the  characteristic  scalelength  of  such  a  spectrum  will  usually  be  larger 
than  that  cf  a  “bump”  spectrum  (which  can  be  at  any  large  wavenumber  kb  =  2:r/4); 
thus,  for  turbulent  spectra  of  the  same  strength,  the  “bump”  spectrum  would  usually 
have  a  greater  effect  on  laser  propagation. 

We  have  also  studied  the  combined  effect  of  a  “bump”  component  superimposed 
on  a  cascade  spectrum.  One  important  conclusion  which  results  from  this  study 
is  that  in  three  dimensions,  even  if  the  amplitude  of  the  “bump”  component  of  the 
spectrum  is  much  smaller  than  the  amplitude  of  the  cascade  component  (for  example, 
by  ten  orders  of  magnitude),  if  the  scalelength  of  the  turbulence  represented  by  the 
"bump”  is  much  smaller  than  the  scalelength  of  the  cascade  turbulence  (say  four 
orders  of  magnitude  smaller),  then  almost  all  of  the  power  in  the  spectrum  is  due  to 
the  “bump”  component.  Furthermore,  in  this  case,  the  effects  of  the  turbulence  on 
the  beam  are  in  effect  determined  by  the  “bump”  component,  and  the  existence  of 
the  cascade  component  (which  may  represent  turbulence  at  very  large  scalelengths) 
can  be  ignored. 
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